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Einstein field equations for anisotropic spheres are solved and exact interior solu- 
tions obtained. This paper extends earlier treatments to include anisotropic models 
^ ' which accommodate a wider variety of physically viable energy densities. Two classes 

b : 

of solutions are possible. The first class contains the limiting case /i oc r ^ for the 
energy density which arises in many astrophysical applications. In the second class 
the singularity at the center of the star is not present in the energy density. The 
models presented in this paper allow for increasing and decreasing profiles in the 
behavior of the energy density. 
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1 Introduction 



Anisotropic spheres in general relativity, in spherically symmetric spacetimes which are 
not static, have continued to attract the attention of researchers since the pioneering work 
of Bowers and Liang—. Some of the interesting recent papers published in this regard 
include the works of Dev and Gleiser«2*i, Herrera et al^, Ivanov^, Mak and Harko^, and 
Sharma and Mukherjee^. The original reason for studying relativistic anisotropic spheres 
was to generate models that permit redshifts higher than the critical redshift Zc of isotropic 
matter. The analysis of Bondi- demonstrated the feasibility and stability of higher red- 
shifts for anisotropic Newtonian and relativistic stars. Anisotropy cannot be neglected in 
stellar clusters and galaxies, in addition to individual stars, as pointed out by Binney and 
Tremaine^, Cuddeford^^ and Michie^. Anisotropy may be an intrinsic feature of boson 
stars as suggested by Dev and Gleiser— and Mak and Harko>£. Sharma and Mukherjee^ 
considered the theoretical possibility of anisotropy in strange stars, with densities greater 
than neutron stars but less than black holes, which are bound stable stellar bodies. 

In a recent treatment Chaisi and Maharaj— found an exact anisotropic solution to the field 
equations utilizing an algorithm of Maharaj and Maartens^^. This solution has a simple 
form: it can be expressed completely in terms of elementary functions and it has a clear 
physical interpretation. In this analysis we extend the Chaisi and Maharaj— solution to 
include anisotropic models which allow for a wide variety of densities which are physically 
reasonable. We present four classes of solution and discuss their physical features. In 
§2 we present the relevant classes of solutions in the form of tables to ensure clarity and 
to highlight the interconnections between the classes. In §3 the physical features of our 
solutions are briefly discussed, and §4 contains the discussion. 
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2 New classes of solutions 



The Einstein field equations to be integrated are 

e-' = 1 1 

r 

r(r — 2m)v' = PrV^ + 2m (2) 
4 

(/i + pr) v' + 2p'^ = — {pr - p±) (3) 

r 

In the above we have defined m(r) = ^ x^/i(x)dx which is the mass function. The radial 
pressure Pr = p + 23/ ^/S and the tangential pressure p± = p — S/ -\/3 are not equal for 
anisotropic matter. The magnitude 5* provides a measure of anisotropy. 

The field equations were integrated by Chaisi and Maharaj— for the energy density 
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l^=^ + k + £r' (4) 



r 



where j, k, i are constants. It is possible to integrate the field equations for other choices 
of the energy density which are physically acceptable. In this section we present three new 
classes of exact solutions to supplement the class found earlier. 



Case 


/i(r) 


I 


:^ + k + £r'^ 


II 


^{l-a)-{p + l)brP-^ 


III 


1 f a-c+{b-d)r^\ ^d[a+br^)-b[c+dr^) 


r2 a+br^ J ^ {a+br^f 


IV 


1 f a-l+6r2+cr*^ , o b+2cr'^ 
r2 \^ a+br2+cr* J ' ^ (^a+br'^+cr'^ f 



Table 1: Energy density functions 



Table Q comprises a list of forms for energy densities that have been studied. The forms 
of energy density fi in Table ^ were chosen so that the Einstein field equations could be 
fully integrated, and the gravitational potentials and matter variables written in closed 
form. The functions chosen for fi in Tabled have profiles which correspond to physically 
acceptable anisotropic spheres, and their simple forms facilitate the analysis of the gravi- 
tational potentials and the matter variables. To complete the integration of the Einstein 
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field equations we also need to make a choice for the radial pressure pr. Clearly a variety 
of choices for is possible; our choice is made on physical grounds. The following is a list 
of forms for the radial pressure: 



Case 


Pr (r) 


I 


^i^-J-lr'-lr') (l 


i?2 ) 


II 


f (a + brP) (l - ^) 


n 


III 


c \ a+br-' J \ J 


n 


IV 


aC 

a+br'^+cr^ \ R? ) 



Table 2: Radial pressure functions 



The forms of Pr selected in Table |21 all reduce to the expression pr = C (1 — r^/i?^)" 
with appropriate choices for the parameters so that the radial pressure is a monotonically 
decreasing function. Table El contains the gravitational potentials e'^ and e^. Table E] lists 
the corresponding matter variables: the energy density fi, the mass function m, the radial 
pressure pr and the tangential pressure p±, respectively. 

We believe that the families of solutions for Cases I- IV presented in Tables EHU are new so- 
lutions to the field equations, apart from particular special cases, for relativistic anisotropic 
matter. Case I contains the Maharaj and Maartens^^ solution {j = i = 0) and the Gokhroo 
and Mehraii solution (j = 0); I2 is defined in terms of elementary functionsJ^. The exact 
solutions in Cases I-IV are amenable to a physical analysis because they have a simple 
form, and in all cases the gravitational potentials and matter variables are given in terms 
of elementary functions. The parameters j, k, and i in Case I are constants. The quantities 
a, b, c, d and p are also constants in Cases II- IV. Note that the constant C corresponds to 
central radial pressure (C =Pr{0),n > 1) 
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Case 




e 






I 


(l-j-|r2-|r^)-^ 


firi-^ exp|/i^. 2(1- 


-Mn+l) y R^J 




II 


ip>0) 






III 


a+br^ 


a — c , - be — ad 

Br c (c + ar ) 2cd exp < 


aCR^ fi r2 ' 
2c(n+l) l^-^ i?2^ 




IV 


a + + cr'^ 




2 \ J 2 
/?2 y 2 


cr I 



Table 3: Gravitational potentials 



3 Physical features 

The gravitational potentials e'^ are finite for all Cases I-IV at the centre r = and at the 
boundary r = R. The functions are continuous and well behaved in the interior of the 
relativistic star. The gravitational potentials e" for all Cases I-IV are continuous and well 
behaved in the interior and finite at the boundary of the star r = R. However we observe 
from Table El that there is a singularity at the centre r = in general for all Cases I-IV 
in the potential e'^. The singularity in e'^ is not present for specific choices of parameter 
values and may be removed by setting 

j = 0, in Case I, (5) 
a = 1, in Cases II and IV, (6) 
a = c, in Case III. (7) 

The gravitational potentials in Table El have the advantage of having a simpler analytic 
form, and they are written in terms of polynomials, rational and exponential functions. 
Consequently the radial and tangential pressures have a simple analytic representation. 

The radial pressure Pr is continuous and well behaved in the interior. Also > in 
the interval (0, i?), regular at the centre {pr{r = 0) = C), and vanishes at the boundary 
{Pr{r = R) = 0) in all four cases. The tangential pressure p± in the four cases has a 
singularity at the centre, but is otherwise well behaved throughout the interior of the star 
and finite at the boundary. Note that the singularity in p±_ may be avoided by suitable 
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particular choices for parameter values. In general the tangential pressure is not zero at 
the boundary of the star {p±{r = R) ^ 0) and does not vanish, as does the radial pressure 
{Pr{r = R) = 0). It is also important to observe that the magnitude of the stress tensor S 
is a nonzero function in general for all Cases I-IV. Hence this class of solutions is generally 
anisotropic and it is not possible for S to vanish and obtain an isotropic limit. 

The energy density fi for all cases contains the limiting case 

/i cx (8) 

as can be verified directly from Table |3 It has been demonstrated that the energy density 
(jHI) arises in isothermal spheres for isotropic matter for both Newtonian and relativistic 
stars by Saslaw et al^^. There are many other examples where the form (jHJ arises in 
relativistic astrophysics. Misner and Zapolsky— proposed that isotropic solutions (with 
k = = i) model the physical configuration of a relativistic Fermi gas for some particular 
value of the parameter j. Another example, given by Dev and Gleiser^ who indicate that 
for particular values of j 7^ and A; 7^ (£ = 0) the energy density /i is appropriate 
for modelling the relativistic Fermi gas core immersed in a constant density background. 
Consequently our exact solutions in Cases I-IV, all containing the limiting energy density 
(jHI), may be used to model a variety of similar astrophysical situations when the anisotropy 
is not negligible. 

It is also important to observe that our approach allows us to consider energy densities for 
which (jSJ does not hold. For this scenario we must utilize the values for the constants in 
©-((Tj). With the help of (O-© we observe from Table IHthat the singularity at r = is 
removed and /i is a continuous function throughout the interval [0,-R]. Figure Q provides 
an illustration of the behaviour of the energy density /i, when ((HI)-© is true, for particular 
chosen values of the constants. The radial distance is over the interval < r < 1. From 
Figure ^ we observe that /i is a well behaved function in the interior of the star and has 
finite values at the center and the boundary. Cases I-IV admit both possibilities of yu' < 
and /i' > 0; that is anisotropic stars with decreasing or increasing energy densities from 
the centre to the boundary, respectively, can be studied. 
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Case 



m 



Pr {C = PriO),n>l) 



P± 



P1 



X 



1 _ £1 



x{w6F(i-J-r-f'-')'(i-S)" 



2nC 



(i-J-F 



II 



i(l-a-6r^) 
(p>0) 



X 



^ (a + 6r*') 



X 



P; + S(¥ + Kf-|)-^) (l-ij)%^(a + 6rT^ 



2a? 



+ (1 - a - 6rP) (1 - a - 6(p + l)r^')} 



III 



1 / a-c+{b-c()r2 



2 I ^Tb^ y 



Pr -r r2(a+6r2) 



2(ad-bc) 



X 1 



r 



X 



aCr'^ ad— be 



aCr 



1 



(a+6r2 



!("'^-bc) I aC f c±drl\ _ V\ 



r 

i?2 



a-c+(b-d)r^ I 



r{c+dr'^) 

f-\_ xl^" _ anCr'^ ( c+dr^ 
\^ B? ) cB? \a+br'^ 



IV 



J_ / a-l+br^+cr* 
\ a+br'^+cr'^ 

2(6+2cr2) 
{a+br'^+cr'^f 



(■^ a+br'^+cr*^ 



aC 
a+br^+cr'^ 



Pr 



X 



4 1 r'^{a+br^+cr'^) (a+fir^+cr"')^ 



2(6+2cr2) 



a+br^+cr'^ 



r^ 

i?2 



X |aCr (l - + J(a - 1 + fer^ + cr^)| 



b+2cr^ 
(a+6r2+cr4)2 



1 _ £1 



aCnr'^ / i rf_ 



n-1 



Table 4: Matter variables 



4 Discussion 



We have presented four new classes of solution which model anisotropic stars. These 
solutions may be broadly divided into two categories in terms of the behaviour of the 
energy density fi at the centre r = 0: the first category contains the limiting case fi oc r~^, 
and in the second category the singularity at the centre of the star is not present. The 
solutions found permit ^' < and /i' > which allow for decreasing and increasing energy 
densities as we move from the centre to the boundary of the star. Note that the solutions 
presented in this paper have the feature that Pr 7^ p± in general so that the anisotropy 
factor S ^ 0; consequently our solutions do not have an isotropic limit. 

We make three observations relating to the physical reasonableness of our model. Firstly, 
the vanishing of the pressure at the boundary is a consequence of the first and second 
fundamental forms; there is no restriction placed on the tangential pressure which may 
be nonvanishing. This feature is also evident in the solutions of Chaisi and Maharaj— , 
Maharaj and Maartens^'^ and Gokhroo and Mehra^. Secondly, the solutions in this paper 
do not satisfy an equation of state. For a realistic astrophysical matter distribution we 
should impose a barotropic equation of state relating the pressure and energy density; this 
will be pursued in future work. Thirdly, we observe that for realistic models the space 
derivatives of the matter field should be negative. This will have the effect of restricting 
the arbitrary constants in the energy density; for example this would place restrictions on 
the relative magnitudes of j, k, i of Case I. 
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Figure 1: Energy density jj plots with singularities removed 
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